Large-amplitude Q n -Q p collectivity in the neutron-rich oxygen isotope 20 O 
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By means of HFB calculations with independent constraints on axial neutron and proton 
quadrupole moments Q„ and Q p , we investigate the large amplitude isoscalar and isovector deforma- 
tion properties of the neutron-rich isotope 20 O. Using the particle-number and angular-momentum 
projected Generator Coordinate Method, we analyze the collective dynamics in the {(Qn) , (Qp)} 
plane. The parameterization SLy4 of the Skyrme interaction is used for all calculations in connec- 
tion with a density-dependent zero-range pairing interaction. Our results show that already for this 
moderately neutron-rich nucleus the transition moments are modified when independent neutron 
and proton collective dynamics are allowed. 

PACS numbers: 21.60.Jz, 21.10.Gv, 21 .lO.Ky, 27.30.+t 
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I. INTRODUCTION 

The dominance of the neutron-proton attraction is 
one of the main characteristics of the effective nucleon- 
nucleon interaction. It can be traced back to the isospin 
symmetry of the force and to the cooperation of the T = 
and T = 1 channels. As a result, isovector collective 
modes have a much higher excitation energy than their 
isoscalar counterparts. On the other hand, the study of 
nuclei away from the stability line has led to the dis- 
covery of phenomena (halo and molecular nuclei, reorga- 
nization of the shell structure) which challenge some of 
the standard nuclear structure concepts. These new per- 
spectives have been mostly provided by the analysis of 
neutron-rich nuclei which, because of the relative weak- 
ening of Coulomb effects, give the best opportunities to 
investigate dynamics far away from the N — Z line. The 
new phenomena are partly attributed to a decoupling of 
the two nucleonic distributions when N exceeds Z by 
far. Another cause invoked is the weak binding of the 
last occupied neutrons in the absence of a barrier. Aside 
from new features of ground-state properties, also some 
interesting effects, such as pigmy resonances, have been 
observed which are clear indications that the excitation 
of isovector modes is easier in neutron-rich isotopes than 
in stable nuclei. The pigmy modes are associated with 
dipolar excitations involving a restricted number of nu- 
cleons and partly for this reason, carrying a small fraction 
of the sum rule only. They also correspond to collective 
vibrations. On the other hand, most of the low-energy 
collective structure remains dominated by isoscalar large 
amplitude quadrupole dynamics, leading in many regions 
of the nuclear chart to ground state deformations. 

The strong proton shell closure at Z = 8 and the al- 
ready good experimental knowledge of the complete oxy- 
en isotope series from 16 O up to the drip line at 24 O 
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makes Oxygen an attractive element for a first inves- 
tigation (not to mention that the theoretical analysis 
is significantly simplified in light systems). In fact, 
the experimental effort on these isotopes has spurred 
a number of analysis based either on the mean-field 
method^ pi M, M, [S3, S3, H or on the shell 
model [la , \2a . \2a . \2l\ . Mean-field investigations of collec- 
tive dynamics properties have up to now exclusively con- 
cerned excitations about the lowest configuration which 
always correspond to zero mass quadrupole moment. An 
important result from both experimental and theoret- 
ical works is the weakening of the N = 20 shell effect 
(both 26 O and 28 O are not bound 0, |, S, 1) along 
with evidences for a local subshell predicted sometimes 
at N = 14 [H or at N= 16 0. 

For these reasons, we have chosen to investigate in 
a consistent way the effect of large neutron excess on 
both isoscalar and isovector quadrupole modes of the 
oxygen isotopes. This first article presents our method 
and gives an illustration for the nucleus 20 O. In order 
to treat small- and large-amplitude collective motion on 
the same footing, we perform our analysis within the 
framework of the projected generator coordinate method 
(GCM) based on Hartree-Fock-Bogoliubov 

(HFB) wave functions. A consistent simultaneous treat- 
ment of isovector and isoscalar dynamics is attempted by 
working in a GCM collective space in which neutron and 
proton deformations are allowed to vary independently. 



II. DEFINITION OF THE COLLECTIVE SPACE 

Let us briefly summarize the different steps involved 
in the projected configuration mixing calculation (see 
Refs. [Il, [29|, [H for more details). 

In the GCM, the collective dynamics is constrained by 
the choice of a variational space, the so-called collective 
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space. For the present study, this space is generated by 
a non-orthogonal basis of particle-number and angular- 
momentum projected HFB states. 



A. Constrained HFB 

The HFB states are obtained from self-consistent cal- 
culations using either a single constraint on the axial 
mass quadrupole moment, or a double constraint acting 
separately on the axial neutron and proton quadrupole 
moments. The resulting HFB states are denoted \q), 
where q labels the collective variables. In the case of 
a single constraint, q is identified with the expectation 
value (Q). The dynamics will be studied along only one 
collective dimension, the total deformation. The relative 
contributions of (Q n ) and (Q p ) to the total quadrupole 
moment (Qt) are then entirely determined by the CHFB 
minimization process. In the case of double contraint, 
the proton (Q p ) and neutron (Q n ) quadrupole moments 
are taken as independent constraints, and the collective 
variable q becomes two-dimensional, q = \{Q n ), (Q p )}. 

As in our previous studies [28|, [2(| the Lipkin- 
Nogami (LN) prescription is used to avoid a spurious col- 
lapse of pairing correlations in the HFB wave functions 
for all values of the quadrupole moment. The HFB-LN 
equations are solved with the two-basis method described 
in Ref. [3l|. The HFB states are chosen to be time- 
reversal invariant, and the single-particle states are eigen- 
states of parity, z signature, and projection of isospin as 
in Ref. [32| . They are represented on 1/8 of a cubic mesh 
of size 11.6 fm with a distance of 0.8 fm between mesh 
points. 

We use the effective Skyrme interaction SLy4 in the 
particle-hole channel [H| together with a surface-peaked 
density-dependent zero-range force acting in the particle- 
particle (pairing) channel. Its strength is taken equal to 
— 1000 MeV fm 3 in connection with two soft cutoffs at 5 
MeV above and below the Fermi energies as introduced 
in Ref. [H. 

As a matter of fact, the one-body potentials associated 
with the standard axial quadrupole operator 



Q = 2i 2 



V 



(1) 



have a numerically pathological behavior at large dis- 
tances. Indeed, depending on whether the deformation 
is prolate or oblate, the constraining field decreases as 
— r 2 along the z axis or in the perpendicular plane, re- 
spectively. When solving the mean-field equations on 
a mesh in coordinate space, the finite size of the mesh 
cuts the constraint at large distances, but there always 
remain attractive pockets at the boundaries of the box, 
which will cither distort the single-particle states around 
the continuum threshold, or even bind them. 1 As soon 



as these states have a sizable occupation, as it happens 
in neutron-rich systems, this will often lead to spurious 
results. The large box sizes that have to be chosen in 
neutron-rich systems as discussed here even amplify the 
problem; hence, the constraint has to be damped at large 
distances from the nucleus. We achieve this by multi- 
plying the constraining operator with a Fermi function 
as proposed in Ref. [36] . We first define an equiden- 
sity surface on which the nucleon density p is equal to 
one tenth of the maximal density /O max - The cutoff takes 
place smoothly over a width equal to a for radii exceeding 
the surface by a distance 5r c . 

The parameters Sr c and a have to be selected such that 
they do not affect the results over the range of consid- 
ered deformations. The main role is played by the cutoff 
distance Sr c which we have set to 5.5 fm, slightly larger 
than the recommended value of Ref. [36[ , while the width 
parameter was set to a = 0.4 fm. 



B. Particle-number and angular-momentum 
projected GCM 

The second step of our method is the construction of 
the collective basis 



\JMq) = ^Pi m P N P Z \q), 



(2) 



by means of projection operators, where J\f = 
{q\Poo P N P z \q) 1 ^ 2 is a normalization factor. Although 
the particle numbers of the intrinsic and projected states 
might be chosen different, we do not make use of this free- 
dom and use N = 12 and Z = 8 everywhere, and drop the 
particle number indices for the sake of simple notation. 

The projection operator P^ fK restores the rotational 
invariance broken by the intrinsic deformation of the nu- 
cleus and provides a state with the good total angular 
momentum J. As we consider only axial deformations 
with the z axis as symmetry axis, the projection of the 
third angular momentum component K in the intrinsic 
frame is equal to zero by construction. We limit our 
study to values of the total angular momentum J equal 
to zero and two. 

In the next stage, we construct collective wave func- 
tions from a linear combination of projected states 



\JMk)=Y,f£ q \JMq). 



(3) 



The weights f£ of the projected HFB states are deter- 
mined from the minimization of the total energy 



(JMk\H\JMk) 
(JMk\JMk) 



(4) 



particle orbitals on a basis such as an oscillator basis, the effective 
cutoff is provided by the long-range behavior of the basis wave 
1 When the CHFB calculations are performed by expanding single- functions. 
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The index k has been introduced to label the several so- 
lutions of the minimization equation which we will as- 
sign to the GCM ground state and to excited states for 
each angular momentum J. The variation of the energy 
(SfE'l — 0) leads to the Hill- Wheeler-Griffin equation 



q' 



E'kZqq') fkq> - °- 



(5) 



Its solution requires to compute the off-diagonal matrix 
elements of the energy and overlap kernels 



H 



(JMq\H\JMq') 



l J qql = {JMq\JMq'), 



(6) 
(7) 



which are evaluated with the techniques described in [28| . 
The calculation of electromagnetic transition of order L 
requires an additional evaluation of the projected GCM 
kernels associated with the electric and magnetic tensor 
operators fl L > M=0 h 



T ML,o]J' = {jMq\f^ L ^\J'M'q'), 



(8) 



where J, L and J' satisfy the triangle inequalities. 

Two comments are in order. First, aside from restoring 
the symmetries broken at the mean-field level, the pro- 
jection also performs a transformation from the intrinsic 
to the laboratory frame. As a consequence, every state 
| J = M = q) is spherical in the laboratory frame of 
reference irrespective of the intrinsic quadrupole defor- 
mation associated with the collective variable q. Second, 
the coefficients fj! are directly connected to the collec- 
tive wave-function g^ q which provide informations on the 
spreading of this wave function over the collective space 
spanned by the variable(s) q. A set of orthogonal collec- 
tive wave functions gf, is given by [33] 
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where (J J ) 1 ^ 2 represents the kernel whose square convo- 
lution gives 

Numerically speaking, the GCM calculations are quite 
demanding. Given the number of grid points selected 
on the {(Q n ), (Qp)} surface (about 25), the number of 
gauge angles to restore the particle numbers (7 for both 
protons and neutrons) and the number of Euler required 
for an accurate angular-momentum projection (12), each 
calculation corresponds to the mixing of about 15 thou- 
sand iV-body wave functions. An additional difficulty 
is related the non-orthogonality of the GCM basis that 
generates near-zero eigenvalues of the overlap kernel T qq i , 
which might lead to spurious states. These difficulties are 
treated with techniques described in Rcfs. [28l. [29l. [30j] . 



III. RESULTS 

In the left panel of Fig. [TJ we show the CHFB en- 
ergy curve obtained with a single constraint on the mass 




FIG. 1: Particle-number projected energy curve obtained 
with a single constraint on the mass quadrupole moment (left) 
and the corresponding path in the plane of proton and neu- 
tron quadrupole deformation (right). The dotted line in the 
right panel corresponds to the line of equal proton and neu- 
tron deformation /32, n = (h,p- 



quadrupole deformation as a function of the dimension- 
less quadrupole deformation fa t 



02,t = \ T7T- 



5 4^Q t 



16tt 3R 2 A 



(10) 



where the radius constant is given by R = 1.2 A 1 / 3 fm. 
Please note that in light nuclei and at large deformation; 
the values of the fa as defined through Eq. [10] might 
be much larger than the generating multipole deforma- 
tion in a liquid drop model. Please note also that, as 
in our previous works, the mean-field energy curves that 
are presented are already particle-number projected, to 
avoid the ambiguities related to the Lipkin-Nogami cor- 
rection to the mean-field energy. In the right panel of 
Fig. [T] we show the path followed by the single contraint 
results in the two-dimensional plane of proton and neu- 
tron deformations 



02,n — 



5 4tt Qn 
16tt 3i? 2 A 

5 4tt Q. p 
16tt 3R 2 Z' 



(11a) 
(lib) 



This path corresponds to the valley that is obtained in 
a double-contraint calculation. Please note that, as long 
as the radii of the proton and neutron density distribu- 
tions are not too different, the use of fa )fl and fa tP re- 
moves most of the trivial scaling of the neutron and pro- 
ton quadrupole moments associated with different values 
of N and Z. With such coordinates, equal neutron and 
proton deformation of a saturating system where the neu- 
tron density is strictly proportional to the proton density 
would lie on the first bisector, represented by the dotted 
line in the right panel of Fig. [TJ Already for not too 
large deformations, the path diverts substantially from 
the line of equal proton and neutron deformations, indi- 
cating that for both prolate and oblate shapes the neu- 
trons are more deformed than the protons. 
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FIG. 2: Neutron (top) and proton (bottom) single-particle en- 
ergies (in MeV) versus mass quadrupole deformation (bottom 
horizontal scales) and mass quadrupole moment (top horizon- 
tal scale). Solid (dotted) lines correspond to positive (neg- 
ative) parity states. The dashed lines represent the Fermi 
energies of protons and neutrons, respectively. 



The Nilsson diagrams of neutron and proton single- 
particle energies as obtained with SLy4 for 20 O are shown 
in Fig. [2] along the path of Fig. Q] Qualitatively, both are 
very similar. For neutrons, however, the upward shift of 
the Fermi energy leads to a particle-hole (ph) excitation 
at /?2.t ~ 1-8 ((Qt) ~ 300 fm 2 ) in which the v p 1 / 2 level is 
exchanged with the v fy i 2 m = 1/2 level. The sequence of 
crossings at the Fermi surface for the protons corresponds 
to a 2p-2h ((3 « 0.6) and 4p-4h proton excitation ((3 ss 3) 
and is close to that obtained for 16 O in Ref. [351 ]. 

In Fig. [3] are drawn the deformation energy curves cor- 
responding to three different calculations. The number- 
projected HFB energies are represented by a thin solid 
curve. It displays a pattern typical for an oxygen iso- 
tope: a well-defined spherical minimum and a steep 
rise of the energy versus deformation, with a shoulder 
hinted at moderate proton deformation. For the energy 
curve projected on J = 0, the behavior is markedly dif- 
ferent. The energy does not vary much over the range 
-0.5 < P 2 ,t < 0.8, or -100 fm 2 < (Q t ) < 200 fm 2 , 
equivalently. This reflects that the Op-Oh configuration 
with respect to the spherical HFB solution is the dom- 
inant component of the J — wave function and that 
this component is present in the mean-field wave func- 
tion over a rather large range of deformations. The part 
of the curve beyond @2,t — 0.8, or (Q t ) ~ 120 fm 2 , re- 
spectively, corresponds to the transition associated with 
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FIG. 3: Particle-number projected mean-field (solid curve) 
and particle-number and angular-momentum projected de- 
formation energy curves of 20 O as a function of the mass 
quadrupole deformation. 



a 2p-2h proton excitation across the magic gap at Z = 8. 
These features have been studied in detail in Ref. [H[ for 
16 0. In particular, it has been shown to explain success- 
fully the very low excitation energy of the first excited + 
state of this isotope. By contrast, the 2 + energy curve 
displays marked minima. The fact that for almost all 
deformations but those close to sphericity the number- 
projected CHFB curve lies significantly above the + and 
2 + curves shows that components of higher angular mo- 
mentum have a non-negligible amplitude in the intrinsic 
HFB states \q). 

In the rest of this work, we consider only neutron and 
proton deformations such that the HFB excitation en- 
ergy does not exceed 20 MeV. The associated zone in 
the {(Q n ), (Qp)} plane corresponds to configurations in 
which the proton IS1/2+ and 1^3/2- orbitals remain filled. 

In Fig. [4] are shown the J = (top) and J = 2 (bot- 
tom) projected energy surfaces. The J = surface dis- 
plays two minima, in the quadrants where both neutrons 
and protons have either prolate or oblate deformations. 
One notes that the energy rises slowly along the axis 
(Qn) — in contrast with the single constraint CHFB 
path indicated by the dash line. A similar, although less 
marked tendency exists along the (Q p ) = axis. This in- 
dicates that the + component in the intrinsic HFB state 
is mostly unchanged along these axes and also that there 
already exists a relative decoupling of proton and neutron 
deformations in this neutron-rich oxygen isotope. 

Figure [5] displays the collective GCM wave functions 
9kq> ^1' ©' f° r t nc lowest J = and J = 2 states 
obtained when restricting the collective motion to the 
one-dimensional path given by the dotted line in Fig. 2] 
On the one hand, the J = ground-state wave function is 
localized in the vicinity of (Q t ) — 0, as expected from the 
potential energy surface. On the other hand, it displays 
an asymmetry favoring prolate deformations. 
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FIG. 4: Contour lines of the angular-momentum and particle- 
particle number projected deformation energy surface of 20 O 
(in MeV) as a function of the neutron and proton deformation 
/?2,t, r — n, p. The right and upper axes show the correspond- 
ing proton and neutron quadrupole moments Q T , r = n, p. 
The upper panel shows the surface for J = 0, the lower panel 
for J = 2. The solid line in both panels represents the line of 
equal proton and neutron deformation ^2,p = p2,n, the dot- 
ted line the isoscalar path followed when a single constraint 
on the mass quadrupole moment Qt is applied. 



The representation of the two-dimensional wave func- 
tion is complicated by the fact that the discretized GCM 
wave function, Eq. ([3]), is defined without a volume ele- 
ment. As a consequence, the value of the collective wave 
function gj. q at a given point q represents the wave func- 
tion times an unspecified volume element, that depends 
on the distribution of the surrounding points. On the 
one hand, this definition of the GCM wave function al- 
lows to optimize the set of basis states and to remove 
redundant states without any difficulty. On the other 
hand, any attempt to plot a GCM wave function built 
from unevenly distributed points might be very mislead- 
ing as any curve used to interpolate these points to guide 
the eye is, strictly speaking, meaningless. A plot of a one- 
dimensional wave function g£ q , as given in Fig. [5j is of 
course also affected, but in one dimension one can almost 
always achieve a distribution of basis states that still 
allows an intuitive interpretation of the collective wave 
function. There are further issues, as the only proper 
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FIG. 5: GCM wave functions of the + ground state and the 
lowest J — 2 state of 20 O for a collective motion restricted to 
one dimension, plotted versus the mass quadrupole deforma- 
tion ,02, t- 



metric is defined by the overlap kernels, not the generator 
coordinates, which furthermore are different for each an- 
gular momentum. These points, which affect merely the 
presentation and interpretation of a GCM wave function, 
but not its physical content, deserve further attention in 
the future. Anyway, for the multi-dimensional wave func- 
tion of 20 O the points that can be chosen without having 
redundant ones in the basis of projected HFB states have 
to be irregularly distributed indeed. 2 As there is no un- 
ambiguous way to plot a two-dimensional wave function 
for the distribution of basis states we have used here, we 
will not attempt to do so, but restrict the further discus- 
sion to the observables calculated from it. However, we 
can state that when adding states with a different proton- 
to-neutron ratio of the deformation to one-dimensional 
path from the constraint on the mass quadrupole mo- 
ment, the dominating contribution to the lowest J = 
and J = 2 GCM states always remains the state with 
/?2 ~ 0.8 located on the one-dimensional path. 

The 2 + energy surface (bottom part of Fig. 2|) is split 
into two wells in the fully prolate and oblate sectors (up- 
per right and lower left quadrants) . In the prolate sector 
the energy surface is markedly deeper. This is a good 
indication that the collective dynamics will favor a pro- 
late deformation of the first 2 + state. The minimum on 
this surface is shifted to the region where the neutron 
deformation is larger than the proton deformation. This 
agrees with QRPA calculations 2j|, as well as with ex- 
perimental evidences [2(J. The properties of the first 2 + 
level are summarized in Table fl] 

In an attempt to systematically add states that cover 
the whole n -/3 p plane to the one-dimensional calcula- 



2 There are in fact several near-equivalent sets which have a very 
different distribution in the (Q n ), (Qp) plane. 
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deformation space 


E 


B(E2 1) 




(MeV) 


(e 2 fm 4 ) 


{<&>} 


3.1±0.1 


80±10 


{(Qn), (Q P )} 


3.2±0.2 


180±50 


Experiment 


1.7 


28± 2 



TABLE I: Properties of the lowest calculated 2 + state in °0 
for the two choices of the deformation space considered here 
in comparison with experimental data taken from Ref. 

tion depicted in Fig. [5j it turns out that the GCM ba- 
sis becomes quickly highly redundant, which is reflected 
by the presence of several small eigenvalues of the norm 
kernel. To avoid spurious states which would lead to 
large numerical inaccuracies on energies and transition 
probabilities, one must avoid a very dense and regular 
two-dimensional mesh. This problem is particularly pro- 
nounced in a light nucleus as 20 O, in which a SCMF 
wave function changes very slowly with deformation, but 
becomes more marginal for nuclei beyond the sd shell. 
Therefore, one has to make several calculations with dif- 
ferent choices for the set of SCMF states, and use the 
variational principle to select the mesh giving the low- 
est energies without giving rise to spurious contributions. 
The uncertainty given in Table [J is a measure of the fluc- 
tuations of the results when using a different set of points 
in the {(Q n ), (Qp)} space before entering a zone of nu- 
merical instability. It is noteworthy that the B{E2) value 
(and its equivalent for the neutron density distribution) 
fluctuate much more in dependence of the chosen set of 
SCMF states than the total energies of the first + and 
2 + states. 

The energy gain from the full isospin dynamics is of 
about 800 keV for both the ground state and the lowest 
2 + state, leaving the excitation of the 2 + state nearly 
unchanged. The experimental value is overestimated by 
1.4 MeV. This discrepancy is in the range of values ob- 
served in a global study of the properties of the first 2+ 
states in even-even nuclei [38| , performed with a method 
similar to the one used. The possible origins of this dis- 
crepancy are also discussed in Ref. 38], and for 20 O, 
they are presumably twofold. The variational calcula- 
tion would be improved for the lowest 2 + state if it was 
constructed from a J z = 2 cranked constrained set of 
wave functions rather than from J z = set as is done 
here. There is also a lack of broken-pair two-quasiparticle 
components breaking time reversal invariance in the vari- 
ational space, which should play a role in this semi-magic 
nucleus, see the discussion of chains of heavier semi-magic 
nuclei in Ref. (38j . 

In d isag reement with the global trend found by Sabbey 
et al. [381 ] . the B(E2) value is in much worse agreement 
with the data than the excitation energy. The effect of a 



larger variational space from a double constraint is even 
going in the wrong direction, leading to an overestima- 
tion of the B(E2) value by a factor 6. The value of 
the ratio between neutron and proton transition matrix 
elements is M n /M p « 2.1 > N/Z. This is less than 
the experimental or the QRPA values calculated with 
the same force [21;]. It indicates that the protons fol- 
low too closely the neutrons when deforming the nucleus, 
which leads to the too large B(E2) value that we obtain. 
Since the QRPA value for the B(E2) obtained with the 
same effective interaction is much larger, 3 this probably 
points out to a deficiency of our variational space rather 
than to a problem due to the effective interaction. The 
presence of broken-pair two-quasiparticle states breaking 
time-reversal invariance in the 2 + wave function could 
reduce the B(E2) value and go into the right direction, 
but the present stage of our method does not allow to 
test the magnitude of this effect. 



IV. CONCLUSIONS 

We have started an investigation of the large- 
amplitude quadrupole dynamics in neutron-rich nuclei 
taking into account the full isospin space. Among our 
initial motivation was the search for stable configurations 
in which neutron and proton display different static de- 
formation. Our results indicate that it is indeed the case 
when the usual dimensionless deformations /32,n and /?2, P 
are used to remove the trivial scaling of the quadrupole 
moments with neutron and proton number. 

We have also looked for signs of a dynamical isospin 
instability which would show up in specific spectroscopic 
properties. At a qualitative level, we observe the impor- 
tance of allowing independent proton and neutron de- 
formation: although the main features of the collective 
wave functions are still the same as those determined by 
the one-dimensional dynamics, the full wave function dis- 
plays a significant extension on both sides of this path. 
Similarly, there is a clear influence on the transition mo- 
ments. Unfortunately, the agreement with the exper- 
imental data is not satisfactory for this light nucleus. 
Our model would probably be improved by enlarging the 
variational space for the 2 + state with the inclusion of 
states breaking time-reversal invariance. The develop- 
ments which will allow this inclusion arc still underway. 
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3 It has to be noted that the QRPA calculations of Ref. [2lll make 
simplifying approximations for the residual interaction, which 
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